PROCEEDINGS 

°F SCIENCE 



o 

O 
(N 
> 

o 
O 

(N 



> 

o 
o 



o 



X 



Wilson Loop in Classical Lattice Gauge Theory and 
the Thermal Width of Heavy Quarkonium 



M. Laine 

Faculty of Physics, University of Bielefeld, D-33501 Bielefeld, Germany 
E-mail: 



laineSphysik . uni-bielef eld . de 



O. Philipsen and M. Tassler* 

Institute for Theoretical Physics, University of Milnster, D-48149 Milnster, Germany 



E-mail: ophil@uni-inuenster . de marcus . tasslerguni-muenster . de 



Abstract: We present an estimate for the imaginary part of the recently introduced finite- 
temperature real-time static potential. It can be extracted from the time evolution of the Wil- 
son loop in classical lattice gauge theory. The real-time static potential determines, through a 
Schrodinger-type equation and a subsequent Fourier-transform of its solution, the spectral func- 
tion of heavy quarkonium in finite-temperature QCD. We also compare the results of the classical 
simulations with those of Hard Thermal Loop improved simulations, as well as with analytic 
expectations based on resummed perturbation theory. 



The XXV International Symposium on Lattice Field Theory 
July 30-4 August 2007 
Regensburg, Germany 



* Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http://pos.sissa.it/ 



Wilson Loop in Classical Lattice Gauge Theory and the Thermal Width of Quarkonium 



M. Tassler 



1. Introduction 

Recently, we made an attempt to properly define a static potential in finite-temperature QCD, 
in the sense of obtaining an object which has a direct connection to the spectral function of the 
heavy quarkonium system The spectral function is related to a mesonic correlator obeying 
a Schrodinger equation in real (Minkowski) time, and the corresponding potential was therefore 
introduced as the real-time static potential. Furthermore, employing resummed perturbation theory, 
the real-time static potential was shown to develop an imaginary part, which induces a thermal 
width for the tip of the quarkonium peak observed in the spectral function [^]. In a subsequent 
work [^], we investigated the extent to which there might be non-perturbative corrections to the 
imaginary part, utilising classical real-time lattice techniques. The purpose of the current note is to 
review the results of ref. |j3|], and also to elaborate on our Hard Thermal Loop improved simulations 
in some more detail than in ref. [|3|]. 



2. Real-time static potential 



The heavy quarkonium spectral function in the vector channel, p{(o), can be obtained using 
the relation 



where C>(f,0) is the mesonic correlator 



C>(?,r) 



dte"''C>{t,0) , 



d'x{\j/(t,x+-)Y^W\i/[t,x 



l^(0,0)7^1^)r(0,0) 



(2.1) 



(2.2) 



Here a point-splitting has been introduced to facilitate a perturbative treatment, and W denotes a 
Wilson line connecting the adjacent operators along a straight path. The dilepton production rate 
from (7^-annihilation in a quark-gluon plasma is proportional to the thus defined spectral function. 

Focusing on infinitely heavy quarks, the correlator can be obtained, up to normalization and a 
trivial phase factor, from the analytic continuation of a euclidean Wilson loop 



C>(f,r)ocC£(iY,r), 



1 



CE{r,r) = — Tr (lV(0,r;T,r)W(T,r;T,0)lV(T,0;0,0)W(0,0;0,r)) 

At f 7^ we can write the time evolution in the form of a Schrodinger equation, 

[idt-V>it,r)]C>{t,r)=0, r=\r\, 



(2.3) 



(2.4) 



which defines the object Vy we refer to as the real-time static potential. 

The simplest estimate for y> comes from perturbation theory. An analytic computation with 
proper account taken of HTL-resummation yields the following result in the large-time limit ^ : 



V>(oo,r) 
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exp(— mnr) 
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(2.5) 
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Figure 1 : The resununed perturbative quarkonium contribution (M denotes the heavy quark pole mass) to 
the spectral function of the electromagnetic current, in the non-relativistic regime [co — 2M) /M <g; 1 [b||. 



The real part corresponds to the standard Debye-screened potential of a static quark-antiquark pair 
at finite temperature. The Debye mass is denoted by /Mq. The imaginary part of the potential 



controls the damping of the correlator C>(f,r), which obeys the Schrodinger equation (2.4). 

The spectral function can now be obtained by inserting the static potential into the Schrodinger 
equation, eq. (|2.4{), supplemented by the usual mass term and spatial derivatives, and employing 



subsequently eq. (2.1). The result is shown in fig. 1. The imaginary part of the potential, encoding 
the Landau damping of the off-shell gluons binding the two heavy quarks together, introduces a 
thermal width to the tip of the quarkonium peak. 

As a next step, we would like to estimate Vy{o°,r) beyond perturbation theory. (In principle, 
p{(o) could be extracted from lattice Monte Carlo simulations by means of maximum-entropy 
and related methods; in practice, this involves many subtleties and, possibly, unknown systematic 
errors. For the current status see, e.g., refs. iQ].) A non-perturbative calculation of V>(oo,r) is 
complicated by the fact that a direct analytic continuation from numerical data for C£:(T,r) is not 
feasible. It turns out, however, that the imaginary part of Vy is formally classical [[I]], and can hence 
be probed non-perturbatively with classical lattice gauge theory simulations, of the type originally 
introduced by Grigoriev and Rubakov 



3. Classical lattice gauge theory simulations 

We start our discussion of the real-time lattice techniques by introducing the framework for 
classical lattice gauge theory simulations [^, which is quite similar to the Kogut-Susskind Hamil- 
tonian approach : 

• The fields are discretized using a 3-dimensional spatial lattice. The time coordinate remains 
continuous. 

• Besides the spatial links Ui, corresponding to the discretized colour-magnetic fields, an elec- 
tric field Ei is defined via the relation Ui{x) = iEi{x)Uj{x), where x = (f,x) and U = dU /dt. 
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• A temporal gauge is chosen. The space of physical states is constrained to gauge field con- 
figurations satisfying the discretized Gauss law, 

G(x) = £ [Ei{x) - P-i{x)Ei{x - i)] - fix) = , (3.1) 

i 

with j denoting a possible colour current, and Pi the adjoint parallel transporter, (x + i) = 

Ui{x)^{x + i)UJ{x). 

The classical approximation for Yang-Mills fields at finite temperature follows by supple- 
menting the phase space just introduced with a canonical time evolution and an average over initial 
conditions with a thermal weight. The weight corresponds to the one in the classical partition 
function. 



Z = [ mimi5{G)e-^'' , H = [lReTr(l - Uij) + \tt{E^' 



(3.2) 



where Uij is the plaquette. The classical equations of motion for the discretized system can be 
obtained by invoking the Hamiltonian principle 5S = 0, and read [Tr (TT'^) = S"'^ /2] [^]: 

Uiix) = iEi{x)Ui{x) , Ei = Y^EfT" , E^{x) = -2ImTr [r« £ Uij{x)] . (3.3) 

A more thorough treatment of the long-range dynamics of hot QCD is possible using the 
so-called Hard Thermal Loop (HTL) effective theory [||], which is obtained by integrating out the 
"hard modes" (with momenta of the order of the temperature) from the system, in order to construct 
an effective theory for the soft modes. To keep the effective theory local, certain on-shell particle 
degrees of freedom need, however, to be added to the effective Hamiltonian [^. Once this system 
is discretized and the classical limit is taken, the properties of the hard modes change, and the 
associated matching coefficient, denoted by m^, needs to be tuned correspondingly [10]. In the 
following we denote the new on-shell particle modes by W{x,v). In a numerical implementation 
the following changes are introduced with respect to the classical setup: 



I. The Hamiltonian obtains an additional part, 

(3.4) 



5H = —y 



4n 2 



where W = T"W" (x, v) describes the charge density of the on-shell modes at x moving in the 
direction v = (1, v). 

2. The velocities v need to be discretised. This can be done, for instance, with spherical har- 
monics [11] or with platonic solids [p^]. Choosing the latter approach, we can replace 
/ dQ.y/4nf{v) —>■ l/NpY,n'^if{v„), where Np is the number of vertices of the polyhedron 
used. The equation of motion of the gauge fields then acquires the source term 

j{x) = {am^f— £ v,Wn{x) , Wn{x) =W{x,Vn) . (3.5) 
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Figure 2: Measurement of the correlator Qi(f,r) [3]. The corresponding potential is shown in fig. 3. 
3. Finally, the new fields also evolve in time, according to the following equation of motion: 

Wn{x)=v', {Ei{x) - i [P,Wn{x + i)-P-iWn{x-l)\^ , (3.6) 

where Ei{x) = [Ei{x) +P^iEj{x — 1)]/2. 
For the purely classical simulations, the required set of initial configurations distributed ac- 



cording to the statistical weight in eq. (3.2) and respecting the Gauss constraint was created using 
the following algorithm: 

1. Pre-generate the spatial gauge links Ui with a Monte Carlo simulation of the dimensionally 



reduced effective theory [13] 



2. Generate the electric fields from a gaussian distribution [cf. eq. ( |3.2| )]. 

3. Project onto the space of physical configurations, satisfying the Gauss law [pi]]. 

4. Evolve the fields using the EOM, and repeat from step 2, until the fields have thermalized. 
In the HTL-improved case there are minor changes, but the essence of the procedure is the same. 

4. The imaginary part of the real-time static potential from Wilson loop dynamics 

To obtain the imaginary part of the real-time static potential, a rectangular Wilson loop of 
spatial extent r = |r| and temporal extent t was measured using classical or HTL-improved simu- 
lations. The measured average over a statistical ensemble of initial configurations, as well as over 
lattice sites and loop orientations, is denoted by Cci{t,r) (a typical result is shown in fig. 2). The 



real-time static potential can then be calculated from eq. (|2.4|), 



Vci(^r) ^ , CUt,r) = i-Tr {wHt)WM) , 



(4.1) 



with Wr{t) denoting a spatial Wilson line of length r. Timelike Wilson lines have disappeared due 
to the use of temporal gauge. The result for Vc\{t,r) is purely imaginary, and is shown in fig. 3. 
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Figure 3: The imaginary part of the real-time static potential from the classical simulation (left panel) 
and from resummed perturbation theory (right panel) |^]. 







N 


anio 


confs 


r = \a 


r = 2a 


r = 3a 


r = 4a 


Simulation 


16.0 


12 


0.0 


200 


-0.060(2) 


-0.156(8) 


-0.246(26) 


-0.319(56) 




16.0 


16 


0.0 


160 


-0.059(2) 


-0.155(8) 


-0.245(22) 


-0.326(48) 




16.0 


12 


0.211 


200 


-0.059(2) 


-0.147(7) 


-0.229(23) 


-0.297(51) 




16.0 


12 


0.350 


182 


-0.030(2) 


-0.064(5) 


-0.096(12) 


-0.118(21) 




13.5 


12 


0.250 


142 


-0.071(2) 


-0.174(10) 


-0.270(33) 


-0.341(97) 


Analytic 


16.0 


oo 


0.0 




-0.0601 


-0.1145 


-0.1507 


-0.1737 



Table 1: Overview of the results in the large-time limit [^1. The results from the classical and HTL- 
improved simulations agree within error bars for amo < 0.25 (at j3 = 16). 



As seen in fig. 3, the predictions from resummed perturbation tlieory and from the classical 
numerical simulations are remarkably similar. At the same time, some amplification of the imagi- 
nary part through the inclusion of non-perturbative (and higher-order perturbative) effects is visible 
in the simulation. The difference between the two results becomes more pronounced at later times. 
In particular, in the large-time limit, a difference between the perturbative and the numerical results 
of up to ^ 100% can be observed (at j8 = 16), cf. table 1. 



5. Conclusions 

The results from the real-time lattice simulations confirm the existence of an imaginary part 
in the real-time static potential, indicated already by leading-order Hard Thermal Loop resummed 
perturbation theory. In fact, non-perturbative and higher order perturbative corrections amplify the 
imaginary part, by up to ~ 100%. The amplified imaginary part widens (and lowers) the quarko- 
nium peak in fig. 1, although the qualitative structure remains unchanged. As a side remark, we 
note that the existence of an imaginary part also leads to strong damping in the solution of the 



Schrodinger equation in eq. (2.4), thus significantly facilitating the numerical determination of the 



spectral function through eq. (2.1 ). 
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